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Networks with a given degree distribution may be very resilient to one type of failure or attack 
but not to another. The goal of this work is to determine network design guidelines which maximize 
the robustness of networks to both random failure and intentional attack while keeping the cost of 
the network (which we take to be the average number of links per node) constant. We find optimal 
parameters for: (i) scale free networks having degree distributions with a single power-law regime, 
(ii) networks having degree distributions with two power-law regimes, and (iii) networks described 
by degree distributions containing two peaks. Of these various kinds of distributions we find that 
the optimal network design is one in which all but one of the nodes have the same degree, k\ (close 
to the average number of links per node), and one node is of very large degree, &2 ~ jV 2//3 , where 
N is the number of nodes in the network. 

PACS numbers: 89.20.Hh, 02.50.Cw, 64.60.Ak 



I. INTRODUCTION 
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Recently, there has been much interest in the resilience 
of real- world networks to failure of nodes or to intentional 
attacks Q, S S S S @. The goal of this work is to 
determine network design guidelines which maximize the 
robustness of the networks to both random failures of 
nodes and attacks targeted on the highest degree nodes 



Networks with a given degree distribution may be very 
resilient to one type of attack but not to another. Con- 
sider the simple seven node network example shown in 
Fig. This network is relatively robust with respect 
to a random failure - only a failure of the central node 
will cause the network to fragment. Thus the probability 
that a random failure will cause the network to fragment 
is only 1/7. On the other hand the network is extremely 
vulnerable to a targeted attack - an attack in which the 
most highly connected nodes are removed first. In this 
simple example the probability that a targeted attack 
which removes one node will fragment the network is 1! 

As shown in Fig. ^> we can modify the network to 
make it more resilient to targeted attack by adding more 
links between the nodes on the periphery of the network. 
With this modification, neither a single node random fail- 
ure nor a targeted attack which removes only one node 
can fragment the network. This increased robustness, 
however, comes with a cost. If we define the "cost" to 
construct and maintain a network with a given number 
of nodes as being proportional to the average number of 
links (k) per node in the network, we see that the cost 
of the original network is 12/7 while the cost of modified 
network is 24/7. So for the additional robustness we pay 
a factor of 2 in cost. 

Our goal then becomes how to maximize the robust- 
ness of a network of size N nodes to both random failures 
and targeted attacks with the constraint that the cost 
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FIG. 1: (a) Example of network with low tolerance to targeted 
attack, (b) Example of a network with much higher tolerance 
to targeted attack but with double the cost. 



remains constant. That is, the number of links remains 
constant but the nodes are connected in a different and 
more optimal way. 

Many real world computer, social, biological and other 
types of networks have been found to be scale free, i.e., 
they exhibit degree distributions of the form P(k) ~ fc _A 
BUmmillEElHGll- For large scale free 
networks with exponent A less than 3, it has been found 
that, if nodes fail randomly, essentially all nodes must 
fail for the network to become disconnected 0, 01 • O n 
the other hand, because the scale free distribution has 
a long power-law tail (i.e. hubs with large degree), the 
networks are vulnerable with respect to targeted attack. 
This raises two questions that we address in this work: 
(i) How can we optimize scale free networks to both ran- 
dom failure and targeted attack and (ii) Are there other 
network types that can be better optimized than scale 
free networks. To this end, we begin our analysis with 
scale free networks and then consider networks with other 
types of distributions. 
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II. OPTIMIZATION METRIC 

The threshold for random removal of nodes for any 
degree distribution, P(k), is Q 
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where kq = (k^)/(k). 

Reference [5| describes how to calculate /* arg , the 
threshold under intentional attack. 

A metric we can use to measure the robustness of the 
network to both random and targeted attack is the sum 
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This is only one of a number of possible metrics we could 
use, e.g., we could have used the product f™ nd ■ /* ar s. 
Our results arc, in general, not dependent on the metric 
chosen. 

Our goal can now be stated as follows: for a network 
of a given number of nodes N, how do we maximize /* ot 
while keeping the number of links constant? 

We can estimate an upper bound for /* ot . We first 
note that the maximum value of /^ and is essentially 1 
which is the case when a small number of nodes have a 
very large degree distribution - as in scale free networks 
with A < 3 or in the simplest case where one node is 
linked to all other nodes. In these cases the probability 
of these critical nodes randomly failing approaches zero 
and the threshold is close to 1. The maximum value of 
jtarg j g obtained in the situation in which all the nodes 
have the same degree, (fc), in which case the targeted 
attack becomes equivalent to random failure and we can 
use Eq.JTJ to find / c targ = 1 - l/((fc) - 1). Our upper 
bound /* ot is therefore given by 
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III. POWER LAW DEGREE DISTRIBUTION 

We first study how to optimize a scale free network 
with a single power law regime by varying the exponent 
A and keeping (k) constant |17|. 

In Fig.H we plot the values of / c rand , / c tar «, and / c tot , 
for a range of the exponent A for a network with N = 10 6 
nodes and (fc) = 3 [18J. For this choice of (fc), the upper 
bound of / c tot is given by /* ot « 1.5 [see Eq.(|3J)]. We find 
that: As A increases, /' arg increases but /* and decreases. 
For A w 2.5 /* ot is optimized but the maximum value of 
ytot ^ 1.04) is small relative to the theoretical maximum 
« 1.5. It is interesting that the network is optimized with 
a value of A about 2.5 which is consistent with the range 
of exponents for many real networks 0, 0, 0, 0, 0, fl3 . 

ElliHEJ. 
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FIG. 2: Random, targeted and total critical percolation 
thresholds for scale free networks as a function of the ex- 
ponent A. 



IV. DEGREE DISTRIBUTIONS FORMED BY 
TWO POWER LAWS 

We next analyze a slightly more complex form for P(k). 
Keeping (k) constant, we consider degree distributions 
which consist of 2 segments each of which is a power law. 
The inflection point at which the distribution changes 
slope we denote by a. The hypothesis is that the first 
power law segment (for k < a) with exponent a will 
contribute to the robustness against targeted attack and 
the second segment (for k > a) with exponent A will 
contribute to the robustness against random failures. We 
determine the relative weights of the two segments such 
that /c-totai is maximized. To maintain constant (k) as 
we change a we again adjust the minimum, m, of the 
distribution. 
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FIG. 3: Total percolation threshold vs the inflection point a 
for distributions P(k) composed of two scale free segments 
with A = 2.5 (for k > a) and for the slope (for k < a) a = 
3, 4, 5, 10, 20, 30 (from right to left). 

In Fig. |3J we plot the values of /' ot as a function of the 
inflection point a for A = 2.5, (fc) = 3 and for various a. 
We see that /' ot attains a maximum value that increases 
with increasing a. Thus for a given A we can maximize 
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/* ot by choosing appropriate values of a and a. 

We can further increase the maximum value of /* ot by 
changing the value of A. In plots (not shown) of /* ot as 
functions of a, for a = 10 and various values of A, we 
find that as A decreases, the maximum value attained by 
/' ot increases. Thus we can maximize the robustness of a 
network with respect to both random failure and targeted 
attack by replacing the original degree distribution by 
one with the same (k) but with two power law segments 
characterized by exponents a and A with a large and A 
close to one (the lowest value of A which yields physical 
results) . 

In these distributions with large values of a the total 
probability in the tail of the distribution is a small frac- 
tion of the total probability, so that there is only on the 
order of one node in the tail and, due to the large value 
of a, most of the nodes have almost the same number of 
links - very close to the minimum m. 



V. DEGREE DISTRIBUTIONS FORMED BY 
AN EXPONENTIAL AND A POWER LAW 

With the insight that the larger the exponent a the 
better the optimization, we now consider distributions 
with the initial power law segment (k < a) of the distri- 
bution replaced by an exponential distribution P(k) ~ 
e _/3fe . As expected we find that for a given (3, at some 
value of a, /* ot is optimized and that the optimization 
increases as (3 increases. 



VI. DEGREE DISTRIBUTION FORMED BY 
TWO GAUSSIANS 

Considering the previous cases, it appears that the op- 
timization strategy does not depend on the fact that the 
initial segment of the distribution is a power law or expo- 
nential. Given that the total probability of the nodes in 
the second segment (the tail of the distribution) is very 
small (of order 1), as discussed above, we now want to 
study the case where the second segment is not a power 
law. We therefore consider here a case where the degree 
distribution consists of two Gaussian segments. 

One Gaussian has its center at k\ and width w\ and the 
second Gaussian has its center at k 2 > k% and width tv 2 . 
The ratio r represents the fraction of the number of nodes 
in the second Gaussian to the total number of nodes. We 
consider cases in which r and k 2 are the independent 
variables and fci must be a dependent variable in order 
to maintain a fixed value of (k). 

In plots (not shown) of the total threshold /* ot in terms 
of the ratio r for various values of &2, we find that the op- 
timal /* ot increases and the optimal value of r decreases 
as the value of k 2 increases. In addition, we obtain higher 
values of the optimal /' ot for smaller values of u>i . This 
fact indicates that the highest value of /* ot is achieved in 
the limit where this width goes to zero. In this limit the 



lower segment tends toward a simple delta function. This 
observation motivates us to study next the optimization 
of networks consisting of two delta functions . 



VII. DEGREE DISTRIBUTIONS FORMED BY 
TWO DELTA FUNCTIONS 

Next we consider the degree distribution that consists 
of two delta functions: 



P(k) = (1 - r)S(k - kx) + rS(k - k 2 ). 



(4) 



As in the case of two Gaussian segments, we calculate 
the total threshold as a function of r and k 2 for a fixed 
value of (k). We obtain analytical expressions for both 
/ c rand and / c targ as follows. 
Using Eq. (JTJ, 
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{k} 2 - 2r(k)k 2 - (1 - r)(k) + rk\ 
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For the threshold for targeted attack, we must consider 
two cases: 

(i)/* ars > r. In this case, after the targeted attack, the 
only nodes that remain have degree k\ . We find 

1 - r 



(k) - rk 2 



(k) - rk 2 - 2 (1 - r) 
{k} (k)-rk 2 -(l-r) ^ 



(6) 

(ii) /' arg < r. For this case nodes are removed only 
from the higher segment and we find 

, targ = (k) 2 -2r(k)k 2 + rk 2 - 2(1 -r)(k) 
Jc fc 2 (fc 2 -l)(l-r) • 1 1 

With the expressions for the thresholds, Eqs. JSJ), JHJ, 
and J7J), we are able to evaluate the total threshold /' ot . 
We can obtain an expression for the optimal value of k 2 
as a function of r by determining the value of k 2 for which 
/' ot is maximized. Based on our results above, we expect 
the maximum will be obtained for r small. Using Eqs. © 
and Q), we find that for small r the optimal value of k 2 
can be approximated by 
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r- 2 ' 3 = Ar-*' 3 . (8) 
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Using this result and Eq. J2J we find, for small r, 



/c t0t = ft 



(9) 



Thus /* ot approaches the theoretical maximum value 
when r approaches, but is not, zero. For a network of 
N nodes, the maximum value of /* ot is obtained when 
r=l/N the smallest possible value consistent with there 
being 1 node of degree k 2 . Given this r the equation 
determining the optimal k 2 is 



ki = AN 2/3 . 



(10) 
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Figure0]demonstrates the behavior of the optimal /' ot 
as a function of r and k 2 . We see that the highest values 
of the optimal /' ot are attained as r approaches zero; and 
for a given small value of r, /' ot is optimized for k 2 from 
Eq. (HJ. 
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FIG. 4: Contour plot of total percolation threshold vs r and 
for distribution consisting of two delta functions with (k) = 3. 

The general nature of our results hold for the metric 
defined in Eq. J2J as well as f° r metrics /* ot defined as a 
linear combination of the random and targeted thresholds 

fc = a fc + °fc : I 11 ) 

where a and 6 allow one to specify for a given network 
the importance to be attached to random and targeted 
attack respectively. The only modification to our results 
for these alternative metrics, is that the prefactor A is 
generalized to 

VIII. DISCUSSION AND SUMMARY 

We develop a strategy for optimization of scale free 
and two-peaked networks against both random failures 
and targeted attacks. To our knowledge, this is the first 
study of the robustness of complex networks to multi- 
ple types of failure/attack. We find that the network 
which approaches the theoretical maximum level of opti- 
mization is generated with a degree distribution which is 



non-zero at only two values: k\ and k 2 . This level of op- 
timization is possible because in order to obtain a value 
of J ( rand which is essentially I we have to wire only 1 node 
with a large number of links. The remaining nodes, all 
with the same degree, provide essentially the same high 
degree of resilience to targeted attack as for the case in 
which all nodes have degree (k). Figure compares the 
level of optimization obtained for these optimized net- 
works two-delta-function networks with the level of opti- 
mization obtained for networks with two power laws seg- 
ments and with the theoretical maximum values which 
can be obtained. 
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FIG. 5: Plots of the optimal /' ot vs Ik) for theoretical maxi- 
mum value (solid line), two delta functions (dotted line), and 
distributions consisting of two power law segments (see leg- 
end) 



The optimal network is obtained by connecting k 2 ~ 
AN 2 / 3 nodes to a single node and all of the other 
nodes except the degree k 2 node are of degree k\ ~ 
(k) - A/N 1 ' 3 ~ (k). 

Subjects for further study include (i) an analysis of the 
static and dynamic properties of the optimized two delta 
function networks which we have identified here and (ii) 
the optimization of complex networks under combined 
random failure and targeted attack. Finally we note that 
the origin of the N 2 / 3 appearing in Eq. I|1(J|I may be re- 
lated to the size of the infinite cluster at criticality for 
Erdos-Renyi graphs [13, 123, HJ] . 

We thank L. Braunstein, S. Buldyrev,and S. Sreeni- 
vasan for helpful discussions and ONR and Israel Science 
Foundation for support. 
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